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Abstract 

We calculate the spectrum of blast wave models of gamma-ray burst sources, for various 
assumptions about the magnetic field density and the relativistic particle acceleration 
efficiency. For a range of physically plausible models we find that the radiation efficiency 
is high, and leads to nonthermal spectra with breaks at various energies comparable to 
those observed in the gamma-ray range. Radiation is also predicted at other wavebands, 
in particular at X-ray, optical/UV and GeV/TeV energies. We discuss the spectra as a 
function of duration for three basic types of models, and for cosmological, halo and galactic 
disk distances. We also evaluate the gamma-ray fluences and the spectral characteristics 
for a range of external densities. Impulsive burst models at cosmological distances can 
satisfy the conventional X-ray paucity constraint Sx/S^ ^ few percent over a wide range 
of durations, but galactic models can do so only for bursts shorter than a few seconds, 
unless additional assumptions are made. The emissivity is generally larger for bursts in a 
denser external environment, with the efficiency increasing up to the point where all the 
energy input is radiated away. 

1. Introduction 

In most gamma-ray burst (GRB) scenarios, the very high initial radiation density leads to 
a relativistically expanding fireball. This is unavoidable if the GRB occur at cosmological 
distances (Paczyiiski , 1986, Goodman, 1986), and is also expected in many galactic models. 
For pair-dominated fireballs, an energetic but very short burst is obtained, whereas in the 
more common baryon-dominated fireballs most of the initial thermal energy gets converted 
into kinetic energy before it can be radiated away (Cavallo and Rees, 1978, Paczyhski , 
1990, Shemi and Piran, 1991), and only a weak burst is expected. However the kinetic 
energy of the coasting baryons can be reconverted into radiation when the ejecta are 
decelerated by the external medium (Rees and Meszaros , 1992, Meszaros and Rees, 1993). 
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This solves the total energy problem; it also results in an observable burst with a longer 
timescale since this is then determined by the dynamic timcscale at the deceleration radius. 
For magnetically dominated fireballs the relativistic outfiow is predominantly in the form of 
Poynting flux (e.g. Narayan, Paczyhski and Piran, 1992, Usov, 1992). In all blast 
wave moves ahead of the contact discontinuity, while a reverse compression shock moves 
into the ejecta, randomizing the directed kinetic energy of the fireball. These shocks can 
accelerate particles which, in the presence of magnetic fields (either frozen-in or created by 
turbulent instabilities) lead to a radiation burst whose total energy is comparable to the 
initial energy input and whose time-delayed duration in the observer frame is comparable 
to that observed. 

The radiation in these models depends on the efficiency of particle acceleration 
behind the shocks and on the magnetic field strength. The details of the particle accel- 
eration process in relativistic shocks are uncertain, as are the mechanisms determining 
the magnetic energy density. We therefore explore various possibilities. Our first aim 
is to discover what range of parameters (e.g. expansion bulk Lorentz factors, external 
density) result in gamma-ray production with an efficiency and timescale appropriate for 
the observed bursts. Secondly, we consider the spectrum of the resulting radiation, and 
particularly the extension of this spectrum towards both softer and harder energies. This 
is interesting, among other reasons, because of the well-known 'X-ray paucity constraint' 
on some observed bursts. In our models, the spectra generally consist of two or more 
components, and are therefore not a simple power-law. The relative intensities in the 
different bands depend on the frequencies at which the breaks between these components 
occur. Observations (or even upper limits) in the X-ray band, or in the optical and UV, 
can therefore help to pin down the parameters of the model. 

Our relativistic blast wave models were originally applied to cosmological bursts. 
However, the same phenomenon would occur, on a smaller scale, when relativistic plasma 
released, e.g., by magnetospheric disturbances on a galactic neutron star runs into the in- 
terstellar medium (Begelman, Meszaros and Rees, 1993). Our considerations are therefore 
also relevant for galactic models of gamma-ray burst sources, if one assumes that suffi- 
cient sources with the appropriate energy are available. We stress, however, that the blast 
wave burst model is insensitive to the specific energy input source. It is almost irrelevant 
whether the initial energy source is a merging neutron star binary, a failed supernova, a 
collapse-induced giant magnetic flare, or any other process with the right total energy. In 
this paper we examine the efficiency, duration and spectral properties of both galactic and 
extragalactic GRB blast wave models. 

2. The Shock Burst Scenario 

2.1 Standard Parameters and Basic Picture 

For a flreball moving into an external medium characteristic of a cosmological scenario, 
the typical parameters are next ~ 1 no cm~^. Eg ~ lO^^E^i ergs, r] ~ 10^773, where Eg 
is the initial flreball energy, 77 = Eo/MoC^, Mq is the entrained baryon mass and n^xt 
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is the density of the external medium into which the fireball expands. The input total 

energy Eo may be initially dominated by either a radiation/pair mixture or by magnetic 
fields, depending on the initial conditions of the burst. The corresponding numbers for a 
galactic scenario would be next ^ 10~^n_3 cm~^, Eo ~ 10^^i?4i ergs, rj ~ 10^r]2 (halo) and 
next ~ 1 no cm~^, Eg = lO^^^^sg ergs, rj = 10^?72 (disk). The bulk Lorentz factor initially 
grows as r(r) oc r, and eventually saturates to the value F/ ~ 77, when the internal energy 
per baryon becomes non-relativistic. Acceptable models for impulsive fireballs generally 
require i] ^ fO^ — fO^ (e.g. Meszaros and Rees, f993) in order for the burst energy 
and timescale not to be too low. The observer-frame deceleration radius r^, and the 
corresponding comoving-frame expansion timescale tex and density of the ejecta material, 
are, for these three scenarios, 

Td ~ lO'W^'Eli'e-'/'v-'/' cm ~ 2.2 X lO^^nlJ/^E^f ^-^/^ry"^/^ cm ~ 

4.6 X 10''n-'/'Eli'e-^/'v2^/' cm (2.1.1) 
te. ~ lO'n^'/'Ell'e-'/^'/' s ~ 2.2 X lO^nlJ/^Ej/^^-^/^^"^/^ s ~ 

4.6n-'/'EU'e-'/%'/' s , (2.1.2) 

Hd ^ lO^no'T]'^ cm~^ ~ lO^noryi cm~^ ~ 10~'^n_3?7| cm~^ , (2.1.3) 

where rj is the bulk Lorentz factor at the start of deceleration and 9 is the opening half- 
angle (if not spherical) along which the fireball expands. The observed (lab frame) burst 
duration is rj times shorter than the comoving dynamic time tex, i-e., 

tL ~ 1 n^'^'Ell'e-'/%'/' s ~ 2.2 n-_l/'El{'e-'/^V2'^' s 

~ 4.6 X lO-^n-'/^Eli'd-^/'rj^'/' s . (2.1.4) 

Generally the deceleration of the ejecta occurs after the bulk Lorentz factor has saturated, 
''d ^ ~ V^oj and also after rh ^ if'To (Meszaros , Laguna and Rees, 1993), so the 
lab-frame expansion time tj, ~ f(iC~^r]~'^ ^ ro/c ~ 10~'^rQ s is of the order of seconds, 
and occurs well after the pairs of the ejecta drop out of equilibrium. (The ejecta become 
optically thin to baryonic electrons at a radius smaller than the deceleration radius, when 
a weak and very brief burst of tt ~ Vo/c ~ 10~* s is produced which precedes the shock 
burst by a time of order ~ ti,). 

A strong relativistic blast wave moves ahead of the fireball. This blast wave 
decelerates as it sweeps up external matter; a reverse shock wave starts to move into 
the ejecta, and becomes marginally relativistic when the expanding ejecta reach a radius 
~ Td- New pairs may be formed in these shocks, whose maximum effect on the opacity 
can be evaluated by estimating the comoving compactness parameter at the deceleration 
radius, assuming all the kinetic energy available is converted into radiation above the pair 
formation threshold. This is r± ~ Lar/meC^AR ~ n±aTAR ~ 10'^ El(^nl^^d^/^r]l^^ ~ 

<^ 1, i.e. both shocks are optically thin to pair formation (and also 
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to baryonic electrons). The spectrum of the radiation produced by the shocks is therefore 
virtually assured to be non-thermal, since the shocked electrons are relativistic and emit 
in an optically thin environment. Indeed, anything other than a non-thermal spectrum 
would be surprising under these circumstances. 

The foregoing discussion, which treats the primary energy production as instan- 
taneous, is self-consistent provided that the fireball is created in a time short compared 
with tL- The actual energy release will, of course, really have a finite duration, and the 
fireball's early development will depend on the time-history of the energy production and 
on the behavior of i]. An interesting possibility, especially relevant to 'cosmological' mod- 
els which invoke a coalescing compact binary, is that the primary energy emerges as a 
magnetically-dominated wind lasting for up to a few seconds (e.g. Narayan et al. 1992; 
Usov 1992,1993) . Some gamma rays may come directly from the inner part of the wind 
(just as they can come from the 'thinning' stage of an impulsive fireball). However, just 
as in the impulsive case, the wind will inflate a cavity, piling up external matter behind a 
blast wave, and may thereby convert its energy into gamma rays more efficiently. 

To explore fully the consequences of non-impulsive (and possibly time-dependent) 
energy release, one would need to introduce further parameters. It is premature to do this 
in an elaborate way. It is, however, worthwhile and interesting to consider a simple (and 
physically plausible) illustrative case when the wind itself is unloaded and almost purely 
electromagnetic, so that its speed is not signiflcantly below c. Suppose this wind turns on 
at t = 0, and maintains a steady isotropic energy flux L for a time tw It will inflate a 
spherical cavity of radius r. If the contact discontinuity at the cavity boundary expands 
with a Lorentz factor T{r), the wind pressure, in the frame of the contact discontinuity, 
is proportional to L/r^F^. This would be balanced by the pressure of the shocked ma- 
terial between the blast wave and the shock discontinuity, which scales as Tig^^ r2. The 
expansion Lorentz factor will be that for which these two forces balance, and this implies 

-r^ _l/2 -1/4 

r oc r / n^J . 

When the wind ceases, at a time tyj after it switches on (measured in the ob- 
server frame), the part of the blast wave directed towards the observer will have reached 
a radius r = T^ct^ and rj will have fallen (owing to the 'loading' by swept- up matter) to 
~ {L/Anc^tf^jUexty^^- At that stage, the configuration would be essentially equivalent to 
a fireball where E = Lt^ and where rj has the value defined above. For these values, we 
find (as consistency requires) that the value of given by (2.1.4) is indeed equal to tw 

Although we present our subsequent discussion in terms of impulsively-produced 
pair, baryon or magnetically dominated fireballs, from the above discussion it applies 
equally to the corresponding wind models at the phase t ~ when they are releasing 
most of their energy. The radiation from an impulsively-produced homogeneous fireball 
rises to a peak after a time ~ rdC~^rj~'^ and then fades. The details of the initial rise 
would be slightly different in the wind case as indeed they would for a (more realistic) 
inhomogeneous fireball that could not be characterised by a single value of ry. In this 
paper we are concerned with the intensity and spectrum of the radiation emitted near the 
peak; the key parameters are then ry, and next- 
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2.2 Shock Structure, Radiation Densities and Timescales 

A crucial question is whether the radiative efficiency of the shock(s) is sufficiently high, and 
this depends on the magnetic field. The synchrotron radiative efficiency of the deceleration 
blast wave and of the reverse shock are near unity if turbulent instabilities at the shocks 
lead to magnetic fields not too far below equipartition. Shocks are unstable (e.g. Ryu and 
Vishniac, 1992), which should result in turbulence, and turbulent field growth seems to 
be efficient in observed supernova remnant and radio source jet shocks. Alternatively, the 
ejecta may be magnetically dominated, or at least contain magnetic fields which initially 
contribute a fraction ^ of the dynamical equipartition value, and subsequently remain 
frozen-in, so oc where V is the comoving volume. In this case the deceleration- 

induced reverse shock is expected to have high inverse Compton (IC) radiative efficiency 
(Meszaros , Laguna and Rees, 1993), for moderate values of the electron Lorentz factor. 
Previously, we had assumed that the IC scattering occurs on the same electrons that are 
responsible for the synchrotron photons, an assumption which is relaxed here. In both the 
turbulent generation and the frozen-in scenario, the main uncertainties are the value of 
the magnetic field at the deceleration radius, and the value of the typical Lorentz factor of 
the radiating electrons, the latter being dependent on the type of acceleration mechanism 
at the shock. However, under the minimal assumption that electrons are accelerated to a 
power law distribution, it is possible to derive reasonable lower limits to the synchrotron 
and IC efficiencies under various parametrizations of the field strength and the minimum 
electron Lorentz factor. 

The detailed structure of the shocks will depend on the magnitude of the field and 
on the efficiency for accelerating electrons to relativistic energies. The strength and Mach 
number of the blast wave and the reverse shock should be difi^erent, as well as the seed 
field density, so in principle the field growth and acceleration efficiency could be different in 
both. The importance of the synchrotron and IC mechanisms will depend on the strength 
of the magnetic field behind the blast wave and the reverse shock, and several possible 
emission region configurations can be considered. 

One extreme case is that where the ejecta have a frozen-in magnetic field left over 
from the explosion phase, at which time its energy density made up a fraction ^ of the total 
initial radiation density. In this case the magnetic field strength in the reverse-shocked 
region at the deceleration radius rd is B^f ~ OA^^/'^E^^^^riQ^^rj^ G . On the other hand, 
if the shocks are subject to instabilities, leading to magnetic field growth which reaches 
some fraction A of the equipartition value with the corresponding post-shock energetic 
particles (as in supernova remnants, or radio sources), then the field at the deceleration 
radius is B^e ~ 4 x 10'^ nl^^rj^X^/'^ G. Notice that the equipartition (or sub-equipartition) 
field value is the same (if A is the same) behind both the blast wave and the reverse shock, 
because the particles behind both must be in pressure equilibrium with each other. The 
corresponding magnetic energy densities ( ergs cm~^) are 

r 6.6 X 10^noA?73 ergs cm~'^ , ( turbulent ); 

Ub ~ S _o 4/3 —1/3 A _Q / X (2-2.1) 

1^ 7 X 10 Uq E^-^ ' ^?7| ergs cm , ( frozen-in ). 
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The parameter A could be different in tfie blast wave and reverse shock (in the equipartition 
case), due to the different physical conditions, while the frozen- in case refers only to the 
reverse shock. 

The acceleration mechanism is probably more efficient when the shock is strong 
and relativistic. Because the external medium is relatively cold and devoid of strong fields 
its sound speed is low, and the blast wave is certain to be strong and ultrarelativistic 
since the bulk Lorentz factor F is large. If the acceleration mechanism does not depend 
on the presence of fields stronger than the minimal seed field expected in the interstellar 
medium, an electron power law can be expected to form behind the blast wave. However 
the strength of any synchrotron radiation from these electrons does depend on the field 
strength. Thus synchrotron emission from the blast wave will be generally negligible for 
the fields B ~ 10~^ G typical of normal ISM, but will be important if turbulent field 
growth occurs. IC scattering, however, only requires the presence of relativistic electrons, 
and may be expected to be important under a wide range of circumstances in the blast 
wave. 

The reverse shock which moves into the ejecta encounters gas which is generally 
hotter than the external ISM, even after adiabatic cooling. For a matter-dominated fireball, 
the sound speed is higher and the reverse shock becomes at most marginally relativistic. 
Similarly, for a magnetically dominated fireball the Alfven speed is high. However, in 
either case, the fireball has a large outward Lorentz factor with respect to the contact 
discontinuity, and the reverse shock may be strong, even if not as highly relativistic as the 
blast wave. Particles are plentiful in the matter-loaded fireball, and in the magnetically 
dominated case particles would be expected due to entrainment or due to sweeping up 
of neutrals from the external medium, which could be accelerated by the reverse shock. 
However, because both the field strength and the particle supply arc different as well as 
the shock strength, it would not be surprising if the blast wave and reverse shock differed 
in their ability to accelerate particles and radiate, unless the interface were unstable and 
mixing occurred on short timescales. 

To cover the possibilities discussed above, we consider several basic shock models. 
The simplest is the "frozen-in" (F) model, which assumes the existence of fields in the 
ejecta at some fraction ^ of the (initial) equipartition value, and no turbulent field growth 
in either the blast or reverse shock. Another model is the "turbulent" (T) model, where 
turbulent field growth is assumed to occur in both shocks leading to fields with an energy 
density which is a fraction A of the particle energy in the shock (A may be different for 
each shock). A third model is the "piston" (P) model, where the ejecta provides pressure 
but the reverse shock is assumed to be an inefficient radiator (either due to low fields or 
poor acceleration), while turbulent field growth and acceleration are efficient in the blast 
wave, which contributes all the radiation. 

If a shock randomizes a large fraction of the total bulk kinetic energy, which 
is comparable to the initial burst energy Eo, the comoving synchrotron luminosity will 
be given by Lgy ~ [Eo/ A'K6'^rf)esyt~^ , where e^y is the synchrotron efficiency relative 
to the other mechanisms (see eq. [2.3.4]). Here we used the fact that in the coasting 
expansion the expanding shell consists, in effect, of ~ A.'nO^rf' causally unconnected regions. 
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each subtending an angle ~ r]~^ at the fireball's centre, and each of conioving dimension 
AR ~ Td/i]] the comoving spectral intensity is Ii, ~ {Lsy/uAR'^) {2i''^ / c^)"^ gym ^^c^ at 
the frequency where self-absorption sets in. The electron Lorentz factor to be used here 
is = max[7TO,, 7a6], where 7^ is the minimum Lorentz factor of the electron power 
law produced in the shock, and 70b ~ is the Lorentz factor of electrons 

that would be radiating photons at a synchrotron frequency equal to the self-absorption 
frequency Vab- The comoving self-absorption frequency is therefore 

^ / L.yC^ri^ .1/3 



7.4 X W^ny\y\esyhmY'\ if 7m > lab\ 

1.1 X W^ny\t^'By^{esy/lQ-^?'\ if 7m < lab- 



(2.2.2) 



The value of e^y and B to be used here is different for the frozen-in and turbulent cases 
(e.g. eq.[2.2.1]). The minimum electron Lorentz factor 7^ of the power law distribution 
behind a shock depends on the bulk Lorentz factor F, and it may exceed F by a factor 
up to mp/rrie ~ 10^. Since the reverse shock eventually becomes marginally relativistic, 
(F^. — 1) ~ 1, while the blast wave (forward) shock achieves a saturation F^ ~ 77 (when 

_ -j^ 2/3 

Tj < F^ ~ 3.3 X 10 Tg ), we will parametrize the minimum electron Lorentz factor 

by 

K, (reverse shock); 



^^^''^'^[Kn, (blast wave), ' ^^.2.3) 

where 1 ^ k ^ 10^ ~ nip/me. The lowest energy electrons emit synchrotron photons with 
characteristic energy i'sy,m ~ ^^^Bj"^^^ ^ 10^ Bk?T'^ ; in the case k ~ 10^ k.^, these are often 
above the self-absorption values (2.2.2). 

The comoving synchrotron cooling time behind a shock characterized by Kgy^sy 

will be 

^ 3.3 X 10^ f 3.1 X lO^no 'A-iK-^r-V' (turbulent); 

^ KsyT^yUB ^ \ 3.1 X lO^n-^/^Eli^^-^Kj^^Tjy^^ s, ( frozen-in ), ^' ' ^ 

which, compared to the comoving expansion time (2.1.2), can be short even for k ~ 1 
in the equipartition case, but is usually long even for k, ~ 10"^ in the frozen-in case. 
The synchrotron energy density can be estimated as Ugy ~ ne(4/3)a"TCUB7^tea;, where the 
optically thin free-flight time across the region AR is essentially the comoving expansion 
time (2.1.2), and the electron density ne behind either one of the shocks is rz-e ~ 
10^nQr]'^r~y (since it is lO^noryf in the ejecta or reverse shock, and 10"^?io?73 in the shocked 
external gas in the blast wave, which has F = Fgy — rf). Thus, while Ugy differs in the two 
field cases, 

^ f 2 X 10-in^/'£;5\/V2/3A«2^F,,r7j/' ergs cm-3, (turbulent); 

\ 2.1 X lO-^nle-'^/^^KlyVsyrjf'^ ergs cm-3, ( frozen-in ), ^ ' ' ' 
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the ratio 

^ = 3 X 10-'nrEl/'e-'/'n%T,yrjl/' , (2.2.6) 
Ub 

is the same in both. 

The comoving IC cooling timescale is given by tie — '^'^eC^ / {^'^tCUsyHic^ ic) ■, or 

1.5 X 10«no-'/'£;5Y/'^'/^A-i«,-2^-^r-ir-S3-^/^ s, ( turbulent ); ^ 7) 

1.4 X 10^So-^^2/3e-^«.t«-^r-^r-S3-^=^/^ s, ( frozen-in ). 

This includes the case when one of the shocks (with Kic, Tic) IC-scatters synchrotron 
photons from the same shock or the other shock with Kgy, Tgy (i.e., if synchrotron and 
IC occur whithin the same shock then Kgy — Kic, = Fjc, or synchrotron and IC- 
scattering occur in different shocks then these quantities differ). The IC cooling timescales 
are, especially for the equipartition case, shorter than the comoving expansion time (2.2.2) 
whenever Kgy or Kic are large (e.g. k ~ mp/nie ~ 10^), and even more so in the blast 
wave, where also Tgy ~ Fjc ~ 10^773. The ratio of the inverse Compton to the synchrotron 
time is 

^ = 4.8 X 10^n-^/^E,Y/^^^/^«7>rc^rri^-V3 , (2.2.8) 

for both the equipartition and frozen-in cases; using "cosmological" parameters, this will 
be short for KgyKicTic '^10^. 

In the case of the blast wave, one must consider also the combined IC cooling 
timescale due to scattering the synchrotron photons that arise in the blast wave plus 
those that arise in the reverse shock and travel through the blast wave towards the ob- 
server. (Note that the bulk motion of the radiating material is only mildly relativistic 
with respect to the contact-discontinuity, so it is a reasonable approximation to suppose 
that ~ I the radiation crosses the contact discontinuity). This timescale is given by 
tic,c = 3meC^/[4crrCKjcrjc('Usy,r + Usy^b)]^ where Ugy^rj Ugy^h ^tre the synchrotron photon 
energy densities due to the reverse shock and blast wave respectively. Since the energy 
loss rates are additive, this combined IC timescale can be obtained from the appropriate 
combination of the timescale (2.2.8), 

'^ic,c — ^ic^rb + '^ic\b ' (2.2.9) 

where rb and bb stand for synchrotron/IC due to reverse/blast and blast/blast, respec- 
tively. Depending on the basic shock model being assumed, some of these components 
may be absent, e.g. in the piston model there are no contributions from the reverse shock, 
and no combined IC scattering. 
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2.3 Spectral Components and Fluences 



The spectrum observed will be made up of a combination of one or two synchrotron spectra, 
and one to three IC-scattered spectra (we consider here only the case where higher order IC 
scattering occurs in the Klein- Nishina regime, so its contribution can be neglected relative 
to the others). The synchrotron spectra will be self-absorbed below Vab given by (2.2.2), in 
which case if 7^6 > 7m ~ the synchrotron spectrum from the blast wave peaks at that 
frequency; the spectra have an energy slope Ii, ^ v^°- above that, where Sa = —{p — l)/2 
and p is the electron power law energy index above 7^ defined by A^(7) oc 7"^. However, 
in the case when «:r ^ 1 so 7^^ > 7a6, the negative slope Sa — — (p — l)/2 starts at 
^sy,m ~ lO^-^Tm > ^ab- The typical values for the non-thermal charged particle power 
index achieved in shock acceleration are p ~ 2 — 4, and taking an average value p ~ 3 
leads to a synchrotron slope above the break Sa ^ —1. For a single field value throughout 
the shock emission region, the slope below this break will be +1/3 while the spectrum 
is optically thin, and -1-5/2 below fab- In practice, however, it is unlikely that the field 
is uniform, and based on the analogy of compact radio sources, which show an energy 
slope ~ below the break, one might expect also here an energy slope Sh below the break 
flatter than either -1-1/3 or -|-2.5. In the absence of absorption, wc shall take a fiducial 
value of Sh — +1/3. The corresponding slope of the power per decade spectrum z/Jj^ will 
be Pi, oc z/"* oc where i — (6, a) for frequencies (below, above) the break. This 

corresponds to the "fiducial" values of the power slope ~ +4/3 below the break and 
Q!a ~ above the break, close to the average values 1 and (e.g. Schaefer, et al. , 1992) 
which have in the past been taken as a guideline. However, one must be aware of the fact 
that there is a considerable spread about these values (e.g.. Band, et al , 1993). 

The lab-frame frequency of the synchrotron break or turnover will be given either 
by the comoving self-absorption frequency (2.2.2), or (for the larger Hgy^sy in either shock) 
by 

^sy,m ^ 10^-^7m' both cases blueshifted by Tj — 10*^773 (the bulk Lorentz factor of 
both the reverse shock and the blast wave, in the observer frame). In the latter case, the 
lab-frame synchrotron turnover frequency is 

_ r 4 X Wnl/^X^/'K%T%ri Hz, ( turbulent ); 

\ 4 X 10'nl/'E-,'/V,/'e/'K%T%4 Hz, ( frozen-in ) , ^ ' ' ^ 

The IC-scattered spectrum will show a corresponding IC break or turnover at a 
lab-frame frequency {4/3)Kfj:f^u^yt ~ 1.3 x lO^SK^^Kf^F^^Ff^rya, or 

^ f 5.3 X W^nl/^X^/^K%KlT%Tl4 Hz, ( turbulent ); 

\ 5.3 X 10«nS/^E,Y/^rJ/^e^/^«^,«f,F^,Ff,r7i Hz, ( frozen-in ). ^ " " ^ 

The IC-continuum slope below the turnover will be the same as for the synchrotron spec- 
trum, and that above the turnover will also be the same as long as the electron slope of 
the IC-scattering electrons pic is even slightly steeper than that of the of the synchrotron 
electrons psy, in the case of scattering of reverse-shock photons on blast wave electrons. If 
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the synchrotron and IC electron slopes are the same, the IC slope above the turnover is 
logarithmically flatter than that of the synchrotron spectrum. This logarithmic flattening, 
if present, will be ignored below, since it does not affect substantially the discussion. Al- 
ternatively, if the synchrotron electrons have a slope psy steeper than pjc, the IC-scattered 
photon energy spectrum has a slope of —{pic — l)/2, which we will take for numerical 
examples to be of order —1 (corresponding to a power per decade or fluence per decade 
slope CKa = 0). 

The total spectrum will then show a number of spectral components due to the 
various synchrotron and IC combinations possible between the two shocks, compatible 
with the assumptions made about the field origin. In the simpler case of frozen-in magnetic 
fields in the ejecta, which will affect only the reverse shock (and assuming that turbulent 
magnetic field generation in both shocks is inefficient, A ^ 1), we have three components: 
one synchrotron from the reverse shock, and two IC, one from the reverse and one from 
the blast wave. These may be labeled (sy, r), (ic, rr), (zc, rh). In the case where turbulent 
field generation is efficient enough to achieve a non-negligible fraction of the equipartition 
magnetic field value in both shocks (i.e., a field larger than whatever frozen-in component 
that may be present) there are five spectral components: two synchrotron , one from 
each shock, and three IC. These may be labelled {sy,r), {sy,b), {ic,rr), (ic,bb), {ic,rb). 
Intermediate situations may exist if turbulent field generation is efficient in only one shock 
but not the other, which can be treated in a similar manner. 

Each spectral component will produce a total photon energy fiuence S ( ergs cm~^),| 
which is some fraction ^ 1 of the maximum bolometric fiuence available from the particu- 
lar shock in question. In the most common case when 7/ < = 3.3 X lO^EK^r^ most 
of the initial fireball energy is carried as kinetic energy of motion of the ejecta, which is 
completely re-randomized in the blast wave and the reverse shock. We shall therefore as- 
sume that the maximum bolometric fluence possible So = Sr + Sb is divided about equally 
between the reverse shock and the blast wave, and is 

So = (sr + Sb)iSo/2) = Eo/AttO'^DI = lO-'^E^-.e-^D^i ergs cm-^ , (2.3.3) 

where ~ 1, ~ 1 are parameters describing what fraction of the total bolometric 
fluence is carried by the reverse and blast wave shocks, and -D^ is the luminosity distance, 
Dl = {2c/Ho)[{l + z)-{l + z)^/^] ~ lO^^h-'^ for ^ ~ 1 in an Einstein-de Sitter universe, 
or its appropriate Newtonian scaling for galactic models. 

For each shock we can define the synchrotron and IC radiative efficiencies Csy, eic 
for individual spectral components as 

~ ^'^ e-sy = 1 1 1 ) (2.3.4) 



tj^f, + tsy ~l" tex ij^Q ~l" ^sy ~l" ^ex 



where the energy-loss timescales have to be evaluated for the particular shock j and the 
particular spectral component (combination of shocks j, /c), e.g. in the IC case, one could 
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have ijc,rr) tic,bb or tic^rb- The corresponding photon energy fluences from each shock or 
shock combination due to each mechanism (giving the different spectral component) are 

Sic,jk — ^j{So/2^eic,jk -I ^sy,j — ^j{^o/'^)(^sy,j -i (2.3.5) 

each extending over an energy range that goes significantly below and above the turnover 
(break), to an extent depending on the upper and lower limits of the electron energies. 



2.4 Gamma-ray, X-ray and Optical Fluences. 



We consider now more specifically the fluences in three frequency bands lS.v jv ~ 1 centered 
around 1.2x10^0 Hz (7 -rays near 0.5 MeV), 3x 10^^ Hz (medium/soft X -rays) and 10^^ Hz 
(0/UV), matching approximately the HETE instrument bands (Ricker, 1992). We take 
the synchrotron spectrum to be given by a broken power law in energy flux, ~ zv**, 
where Sj = [sfe, s^] for v <> Vgyt , in terms of the observed lab-frame Vsyt of eq. (2.3.1- 
2.3.2). The IC-scattered spectrum will also be a broken power law with the same energy 
indices, I,^ oc while the corresponding power per decade (and fluence per decade) 
spectra have slopes — vly oc v'^^ where ai = Si -\- 1. In a double- logarithmic plot the 
IC spectrum is shifted with respect to the synchrotron spectrum upwards and to the right 
by a factor (4/3) k^F^. (We assumed, for simplicity, that the scattering electron energy 
index p satisfles ^ {j> — l)/2, although this is not necessary). The 7-ray fluence in the 
~ ~ 0.5 MeV band will then be given by the sum (or in practice by the largest) of the 
various IC and synchrotron spectral contributions in that band. If self-absorrption effects 
can be ignored, this is 

-^7 = XI ^n,jkSn,jk {T^-y/T^ntJkT' ergS Cm~^ (2.4.1) 

n,jk 

where n stands for synchrotron or IC, nt stands for the turnover syt or ict, jk stands for 
the shock combinations rr, bb or rb (where relevant), the constant A^jk ~ 1/few is taken 
for numerical estimates to be ^ = 0.2, and Vntjk stands for the n-mechanism turnover 
frequency of the shock combination jk. Similarly the X-ray fluence in the u ^ — 3 keV 
band is given by 

Sx^^ An,jkSn,jk (I'x/i^ntjk)"' crgs cm~^ (2.4.2) 

n,jk 

and the 0/UV fluence in the ~ z^^ :ii 10^^ Hz band is given by 

Su = ^ An,jkSn,jk {i^u/i^nt,jk)°" ergs cm~^ (2.4.3) 

n,jk 

The above are without self-absorption, and generally we take ctb, eta — 4/3,0, although 
observationally other values are also expected. However, in many cases the spectra are 
more complicated than (2.4.1-2.4.3) because self-absorption becomes important at the 
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lower energies. For Vab < ^m-, the slope is aa for v > 1/^., it is 0:5 (~ 4/3) for Vab < < i^mi 
and 3 for < fab (e.g. 0, 4/3, 3 for decreasing energies). For > ^^m, it is for 
ly > Uab, it is 3.5 for z/^ < z/ < and 3 for z/ < z/^ (e.g. 0, 3.5, 3 for decreasing 
energies). In these cases, an individual synchrotron component has two breaks, and the 
corresponding IC-scattered component will also reflect this fact. Notice that the power in 
a certain band (e.g. eqs. 2.4.1-2.4.3) depends on the slopes a^, but the power at the break 
energies Snjk (eq. 2.3.5) is independent of these slopes. Thus, unless the bands are very 
far from the break energies and the slopes are very different from those assumed here, the 
band fluences should not be too different from those estimated with the fiducial slopes 4/3, 
0. The limiting 7-ray sensitivity of BATSE (e.g. Fishman, 1992) is ~ 10~^ ergs cm~^, 
while for bursts of 10 s that of HETE (Rickcr, 1992) is ~ 3 x 10"^ ergs cm'^, and 
in the X-ray and UV bands the sensitivity limits are ~ 8 x 10~^ ergs cm~^ and 
Su ~ 10~^ ergs cm~^. 

The calculated positions of the spectral breaks, and the slopes of the power-law 
segments of the spectra, are shown in the figures. In these figures, the breaks are depicted 
as being sharp. In reality they will be smoother, not only because an individual electron 
emits broad-band radiation but also because the observed flux comes from different parts 
of the fireball whose doppler shifts have a spread of order 2 (from ~ to almost 2r]) . 

3. Results 

3.1.1 General Features 

In general the fluences scale as D~^, and as a positive power (usually not unity, due to the 
efficiency factors) of the intrinsic total initial energy Eg- They also scale as a negative power 
of the beaming half-angle 6, which is generally not 9~'^ (again due to the 6 dependence 
of the efficiency and turnover factors). The fiuences also scale as a positive power of the 
external density no, generally between 1/3 and 2/3. They also depend on ry; it is the burst 
timescale, however, which is the most sensitive to this parameter. 

In the cosmological case we adopt Eg = lO^^E^i ergs and D = 10^®D28 cm. For 
no = 1 the burst duration for r]s = listL = 5s and the maximum photon bolometric 
fiuence is 5"^ ~ lO~^-E'5i-D~^0~^et, where et is the total radiative efficiency which de- 
pends (usually nonlinearly) on E^IjUq, etc. For smaller rj and/or smaller no the observed 
timescale gets longer, and vice versa. 

In the galactic halo case we use, e.g., E^ = 10^^ ergs, D = 1.5 x 10^^ cm = 50kpc, 
while for the galactic disk we use generally Eg = 10^^ ergs, D = 3 x lO^-*^ cm = 1 kpc. 
The maximum bolometric fluences Sm are (in principle) similar to the cosmological case, 
for this choice. However, especially in the frozen-in field models, which are inefficient 
radiators, higher Eq are occasionally needed to reproduce the observed 7-ray band fiuences. 
Compared to the cosmological case, the durations are comparable for the halo (where 
typically no = 10~^), but for the galactic disk (where typically no = 1) they are shorter 
(see eq.2.1.4). The typical external densities are the interstellar values expected for the 
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halo/disk case, unless the burst occurs inside a denser pocket of material, e.g. a wind 
or nebula that preceded the event. For the lowest interstellar densities, acceptable 7-ray 
fluences require generally large k ~ 10^, while for higher densities no 1 k 40 is 
required. 

The spectrum depends on the value of k, and A in the shocks, and we generally 
assume that k and A are the same (or do not differ much) for both shocks. We shall 
consider values of IQ-^ ^ (A,0 ^ 1, IQ-^ ^ no ^ 10^, k = 1, 40, 10^, 1 ^ r/ ^ 10^ and 
slopes ai = 4/3, in most of the illustrative cases discussed below. 

3.1.2 Frozen-in Field Models (F) 

In this model magnetic fields are assumed to have built up in the ejecta at the beginning 
of the fireball's history to some fraction ^ of the equipartition value with the total initial 
disposable energy E^. Magnetic fields in the blast wave and reverse shock due to turbu- 
lence, etc., are assumed to be weak, so that synchrotron radiation is important only in 
the reverse shock, where the frozen-in field is characterized by ^, but IC scattering in the 
blast wave is important, because it efficiently accelerates electrons to high 7. This scenario 
has only three radiation components. It is also the least efficient, the synchrotron and IC 
radiative efficiencies being generally smaller for ^ = 1 than for the A = 1 shock-turbulent 
or piston models with the same external density (which we discuss below). This is because 
^ = 1 means equipartition at the original explosion, but then the field energy decreases adi- 
abatically with the expansion; even after compression (by a modest factor) in the reverse 
shock, it is well below equipartition with the bulk kinetic energy (see eq. 2.2.1). 

At cosmological distances, if one uses the low value k = 1 in both shocks, and 
9 = = 1, the fiuences in all three wavebands are below the detection threshold for 
D28 = I5 and are of course smaller for ^ < 1. However the situation is better for the case 
when K S> 1, especially if the burst is beamed. For 9 = 10~^, k = 10^ and no = ^ = 1, the 
7-ray fluences are well above threshold, and the major break is in the neighborhood of the 
7-ray band, so that the X-ray fiuence is significantly lower, for a range of burst durations 
comparable to that in observed GRBs. A spectrum for this frozen cosmological (FC) model 
is shown in Fig. la, where the individual components are also indicated. The frozen-in 
models are easily detectable up to D28 ^ 1 for values of ^ not too far below unity. For the 
same parameters but ^ ^ 10~^ the 7-ray fluences would drop below detectability, unless 
-D28 ^ 10~^ or the external medium is denser. For instance, ^ = 10""^, no = 10^, D28 = 1 
yields large 7-ray fluences and reasonable X-ray fluences. 

At galactic halo and galactic disk distances, because of the lower total energy, 
the deceleration radii ra are generally smaller. A burst of given duration requires a smaller 
T] (for a given next)- The magnetic fields are weaker, and the cooling times are longer 
while the comoving expansion time is shorter. Therefore the efficiencies are lower than 
in the cosmological case, and their relative values change. Thus, in the galactic cases 
the synchrotron efficiency generally exceeds the IC efficiency, even for large values of «; ~ 
10^ (unless the burst occurs inside very dense clouds no ^ 10^ cm~^). The spectra are 
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consequently different from those of the same model at cosmological distances. Because 
of the intrinsically lower radiative efficiency of the frozen-in models (compared to the 
turbulent and piston models), detectable fluences are obtained only for relatively large (by 
galactic standards) total burst energies, e.g. Eq — 10^^ ergs and higher for the halo or 
Eo = 10^^ ergs for the disk. Acceptable 7-ray fluences are obtained mainly for beamed, 
high K cases, e.g. 9 = 10~^, k = 10^ and ^ = 1. A spectrum for the frozen halo (FH) case 
is shown in Fig. 2a, for these parameters, uq = 10~'^ and two values of r] corresponding 
to 12 s and 0.5 s durations. A spectrum for the frozen galactic disk (FD) case (no = 1) is 
given in Fig. 3a, for 5 s and 0.25 s. In both halo and disk cases the major break is below 
all three observational bands, so the fluences at UV, X-rays and 7-rays are all comparable, 
Sx/Sf ~ Su/S^ ~ 1. 

The general behavior of the spectra for bursts of various durations is shown, 
for the frozen-in cosmological (FC) model in Fig. 4, for the same parameters as Fig. 
la but variable duration tj,. The three break energies tend to be softer for the longer 
total durations, and also the fluence levels within the same plateau are smaller for longer 
durations. However, for a fixed energy band the fiuence can also increase rather steeply as 
the duration is taken longer, if one of the breaks moves across that energy band. 

3.1.3 Turbulent Field Growth Models (T) 

In this model turbulent magnetic fields are assumed to build up to some fraction A of the 
equipartition value in both the reverse and blast wave shocks (not necessarily the same in 
both, although here we assumed them to be equal). As a result there are generally five 
different spectral components, although often one component dominates over several of the 
bands 7, x, u. For values of k ~ 1 in both shocks, i.e. 7^ ~ F, the synchrotron efficiency 
is generally larger than that for the IC process. For larger k the importance of IC is 
larger, and in the cosmological case (where Eo is larger so is larger and Bd smaller) 
it dominates the radiative efficiency. At galactic distances, Eo and are smaller and 
synchrotron remains mostly dominant, even when IC is important from the spectral point 
of view. 

For cosmological distances, observable 7-ray fiuences are obtained for « = 1 if 
A ^ 10~^, and for k ~ 10^ even for values A well below 10~^. However, for A ^ 10~^ 
this model is also extremely efficient in producing radiation at other frequencies below the 
7-ray band. A turbulent cosmological (TC) spectrum for standard parameters uq = 1, 
K = 10'^, A = 10~^, 9 = 10~^, ?73 = 1, = 5 s is given in Fig. lb, showing the five spectral 
components. It is seen that, even for this low field case (10~^ of the equipartition value) 
the 7-ray fiuence is high. This is true also for higher A values, but the X-ray to 7-ray 
fluence ratio Sx/Sj is then closer to unity, whereas for the lower values such as A = 10~^ 
this ratio is of order ^0.03. It is worth noting that the spectrum for the frozen-in model 
with ^ = 1 plotted in Fig. la is essentially the same as that for the "turbulent" model at 
the same distance with A = 10~^ (in Fig. lb). The reason for this is that for this choice 
of A and ^ the flelds in the shock at deceleration have the same magnitude (eqs. 2.2.1), so 
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the dominant reverse-shock synchrotron , reverse-shock IC and blast wave combined IC 
are the same, while the blast wave synchrotron and blast wave IC spectra are below the 
other three components. 

For galactic halo and disk distances the turbulent model 7-ray fluence is generally 
high, even for densities somewhat lower than the standard interstellar values (no = 10~^ 
halo, no = 1 disk), especially at the higher k values. As in the galactic frozen- in models, 
due to the lower field value the breaks occur at relatively low energies and the flat (slope 
0) spectrum extends to energies lower than in the cosmological case. As a result, there is 
an increased tendency to produce a relatively large X-ray fluence, both the 7- and X-ray 
bands falling above the main break. There are some exceptions to this situation, e.g. in the 
disk case for 6 = 10~^, X = 10~^, k = 10^ the X-ray fluence is well below the 7-ray fluence 
for short bursts r] = 10^, = 0.2 s, but they are comparable for the longer bursts Is, 
and for longer bursts the UV fluence can also be significant. The spectra characteristic of 
such short and long duration halo and galactic disk turbulent models are shown in Figs. 
2b and 3b. 

The spectra of bursts as a function of burst duration is shown for the turbulent 
cosmological model TC in Fig. 5, for fields in the shock which are A = 10~^ below 
equipartition. These spectra show the same break softening and plateau fluence decrease 
with increasing duration as seen in the frozcn-in case, but less pronounced. The reason for 
the difference is the different i] dependence (or II dependence) of the field B at the shock 
(eq. 2.2.1). For this choice of ^ = 1 (frozen) and A = 10~^ (turbulent) the fields (and 
spectra) are the same at tz, = 5 s (r/ = 10^), but differ at other durations (other rj). For 
this low A, the two additional spectral components in the turbulent model just begin to 
become noticeable at the longest durations. The two lower energy breaks are due to the 
reverse shock synchrotron and IC, while the two highest are the blast wave self-IC and 
combined IC (the blast wave synchrotron is just visible around 10^ s near 10^^ Hz). The 
MeV band is in the neighborhood of the reverse IC component (redshift is not included in 
the cosmological plots shown). 

In Fig. 6 we show the turbulent halo model spectra as a function of duration. In 

this case the two main breaks are the reverse and blast synchrotron breaks, the MeV band 
falling near the blast synchrotron break. The turbulent galactic disk model is similar, but 
the MeV band is farther up on the plateau above the blast synchrotron break, while the 
reverse synchrotron is closer in fiuence and energy to the blast synchrotron component. 

3.1.4 Piston Models (P) 

The physics of these models is the same as in the turbulent field growth cases discussed 
before, except that there are no radiative contributions from the reverse shock, only from 
the blast wave. The spectra have only two components, synchrotron and IC. The relative 
importance of the blast synchrotron and IC contributions are changed, relative to the two- 
shock turbulent growth model, because the much more abundant reverse shock photons 
which previously acted as seeds for the blast wave IC are now absent. For this reason, the 
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blast wave synchrotron spectrum dominates in almost all cases. The spectrum being so 

simple, the value of the fluence at each band is determined simply by the slopes assumed 
and by whether the observed band is above or below the synchrotron break. It is therefore 
strongly dependent on the field strength A and the particle minimum energy k. 

The cosmological piston (PC) model has a high bolometric radiative efficiency 
for a wide range of A and k. The highest 7-ray fluences (for a slope above the break) 
are obtained when the break is below the 7-band, which requires A and/or k not to be 
too large. Values of re ~ 10^ are too large, but high 7-fiuences are obtained for k ^ 10^, 
for a range of burst durations. When k is much lower than this, and/or A is very low, the 
X-ray fluence can become comparable to the 7-ray fluence. The spectrum of a cosmological 
piston model and its components is shown in Fig. Ic. 

The piston models in the halo and galactic disk (PH, PD) have a lower radia- 
tive efficiency than the cosmological piston case, due to the lower energy density and 
corresponding lower magnetic field in the blast wave. The total energies required to ob- 
tain acceptable 7-ray fluences are comparable to those in the turbulent growth model, 
Eq = 10'^^ ergs (halo) and = 10"^^ ergs (disk). Due to the lower held the break energy 
appears at low frequencies, typically below the UV band. High bolometric fluences are 
obtained for large k ~ 10^, giving spectra whose fluences in the 7-ray and X-ray bands are 
comparable for a high energy slope of for durations tL ^ 1 s, but ratios Sx/S^ ^ 0.03 
are obtained in the short bursts (ti, ^ 1 s). 

In Fig. 7 we show, for a fixed density, the piston cosmological (PC) spectra 
expected as a function of the burst duration, for A = 10~^. The piston models have 
only blast wave spectral components, and in this case the MeV band falls near the blast 
synchrotron break, while the blast IC component is seen at the highest energies. 

In Fig. 8 the piston disk case (PD) is shown, where again the MeV band falls 
nearest to the blast synchrotron break, while the blast IC component is not sufficiently 
strong to show up. The corresponding piston halo model shows similar properties. 

3.1.5 Gamma-ray and X-ray Fluences 

The spectra previously discussed refer to 'standard' densities, and it is interesting to con- 
sider other values as well. In the density-duration parameter space, we have calculated 
the 7-ray band Sj fluence contour levels and also the ratio S^/Sx, these two bands being 
centered around 0.5 MeV and 1 KeV. These are the most likely to be of immediate use in 
constraining models. Similar contour levels for the S-y/Su ratio can also be made but are 
omitted here for brevity. The S-y contours start at the approximate BATSE threshold of 
10~^ ergs cm~^ and higher values, while the Sj/Sx contours go from 1 to 10^. Of course, 
other values also occur, but these delimit the most interesting range between comparable 
7 to X-ray ratio and the usual "X-ray paucity" value of Sx/Sj ^ few percent. For each 
model these two sets of contours are shown simultaneously in Fig. 9. The left column shows 
the cosmological models (frozen, turbulent and piston from top to bottom), the middle 
column the halo models and the right column the galactic disk models (same distribution). 
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The frozen and turbulent cosmological and S^/S^ contours behave fairly similarly in 
the n,tL plane, not surprisingly since the field values were chosen approximately similar. 
The higher fluences are in the n ^ 1, tL ^ 5 s upper right quadrant, but detectable Sj 
is obtained over much of the range 10~^ ^ tL ^ 10^, 10"*^ ^ n ^ 10^. The region where 
the S^/Sx ratio is larger than 30 is also in the upper right quadrant, so longer bursts 
are better at satisfying the X-ray paucity constraint. The cosmological piston contours 
are much simpler, due to the simpler (essentially one component) structure: burst in a 
medium of n ^ 10^ and with tL ^ 1 s are all above threshold, and those at all densities 
with ti, ^ 12 s satisfy a gamma/X ratio ^ 30. Here it is the short bursts that satisfy 
X-ray paucity. This difference between frozen/turbulent and piston is due to the fact that, 
for this particular choice of A, ^ and k the former have the 7 band halfway up the break 
slope, while the latter has it on the plateau above the break, and the break softens with 
increasing duration. The details of the behavior, it must be stressed, are a function of the 
choice of parameters. 

For the halo models (middle column of Fig. 9) the frozen models are above 7 
threshold everywhere except at long durations and low densities, higher S-y occurring for 
shorter durations and higher densities. For the frozen model the X-ray paucity ratio of 30 is 
violated almost everywhere, while gamma/x ratios ^ 10 are obtained for 5 s, n ^ 10^. 
The turbulent and piston halo are above 7-threshold essentially everywhere, but turbulent 
satisfies a gamma/x ratio ^ 30 for n ^ 10~^ and durations 0.2 s ^ tL ^ 5 s, while the 
piston does it for all densities and tL ^ 5 s. 

For galactic disk models (right column) the frozen models behave similarly to the 
halo case (middle column), but the region where the gamma/x ratio is larger than 10 is 
smaller (no ^ 10^, 10^ ^ n ^ 10^). The turbulent disk has a somewhat larger region in 
the low density, short time quadrant where this ratio ^ 10 is satisfied, and the piston disk 
satisfies it for all densities and ^ 1 s. 

Also shown in Fig. 9 are the lines below which impulsive bursts are not possible, 
shown as full straight lines running diagonally from top left to bottom right in the halo 
and galactic disk cases. In the cosmological cases these lines are outside the figures for 
the range of values used, so this restriction does not apply. Short bursts below and to 
the left of this line would require, in the impulsive burst limit (energy deposition time ^ 
burst duration ti), a value of F = 77 which exceeds the maximum value compatible with 

dynamical requirements, ry > F^ = 3.3 x 10^£'|2"^rg "^^^ (e.g. Meszaros , Laguna and Rees, 
1993). Bursts below this line, therefore, are only consistent with an interpretation where 
the timescale tL is not given by the dynamics as a function of 77, via eq.(2.1.4) (see also 
below) . 
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4. Discussion 



We have calculated the radiation spectra of gamma-ray burst sources arising in the blast 
wave produced when relativistic ejecta encounter the external medium. While there are a 
number of uncertainties concerning the field strength and the particle acceleration efficiency 
in the blast wave and reverse shocks, a simple parametrization has been used to determine 
the range of physical conditions that can be encountered. For models where turbulent 
magnetic field growth occurs in the shock the field strength is characterized by A, which 
is the ratio of the magnetic to the particle energy in the shock, while frozen-in fields are 
characterized by ^, which is the corresponding ratio at the time of the initial impulsive 
event. The particle acceleration is characterized by the ratio k of the minimum electron 
Lorentz factor in units of the bulk Lorentz factor of the burst. The range of parameters 
explored is that which reproduces the observed gamma-ray fluences and the overall burst 
durations, 10~^ ergs cm~^ ;^ ^ 10""^ ergs cm~^ and 10~^ s ^ tL ^ 10^ s. Temporal 
substructure on shorter timescales is expected, since the external medium is likely to be 
lumpy and the shocks are expected to be unstable. These gross features of GRBs are easily 
explained by the impulsive energy deposition model with physically reasonable values of 
these parameters, within the context of the models discussed here. 

A natural explanation for the non-thermal, power-law spectrum seen in the bursts 
in the gamma-ray range is provided by the synchrotron and inverse Compton radiation 
losses expected from the power law relativistic electrons accelerated in the shocks. Also 
the presence of breaks by about one unit in the power law finds a natural explanation in 
this model. These are associated with the minimum energy to which relativistic electrons 
are accelerated. Steeper breaks can also arise from synchrotron self-absorption, which we 
have evaluated here for the case of relativistically expanding sources. This is important in 
some cases, but it appears usually below the UV range. 

The breaks do not appear at a preferred energy, the latter depending on the 
parameter rj and external density as well as k. For a fixed assumed fraction of the equipar- 
tition energy density A, the magnetic field at the shocks depends on rj and the break 
energies depend on B and powers of F. Since for the reverse shock Fj. is always of order 
unity, the blast wave F;, = ry, the rj (or t^) dependence of the break energy is stronger for 
the blast wave components than for the reverse components. The softening of the 'main' 
piston spectral break (due to the blast wave) with decreasing rj (increasing baryon loading, 
or increasing burst duration t^) is more pronounced than the softening of the 'main' break 
(due to the reverse shock) in the turbulent and frozen-in models, but it is comparable to 
the softening of the third, combined IC break (due to the blast wave) in these models. 

Specific constraints are obtained by comparing the fluence in different parts of 
the spectrum. In particular, the usually invoked X-ray paucity constraint S^/S^ ^ 0.03 
(e.g. Laros, et al. , 1984), if applied universally, provides stringent restrictions on the 
type of models which are acceptable. Only the cosmological models are able to satisfy this 
condition over the entire range of durations 120~^ s ^ tL ^ 10^ s: turbulent cosmological 
models are satisfactory for relatively low values of A ^ 10~^, piston cosmological models 
are satisfactory over a wide range of values of A, and frozen-in cosmological models require 
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relatively high values of ^ 10~^. On the other hand, the galactic models produce lower 
magnetic fields in the shocks, and a direct consequence of this is that their breaks occur 
generally at energies below the X-ray band. This holds for impulsive bursts, where is 
determined by 77 < = 3.3 x IO^eK^Tq (eq.2.1.4), for durations ti ^ 1 s. Galactic 
models would therefore not produce 7-ray breaks, and would also predict rather high 
Sx/S'y ~ Lx/L^ ~ 1 fluence ratios violating the X-ray paucity ratio as defined above, 
except for very short bursts with (impulsive) duration ^ ^ given by eq.(2.1.4). A 

possible way around this might be if the bursts are not impulsive, i.e. the original energy 
input has an intrinsic, longer timescale which is identified as the burst duration, e.g. it is 
tyj for a wind model, or taccr for an accretion model. In Fig. 9, impulsive models can only 
exist above and to the right of the diagonal solid line (in the cosmological models, this 
line is outside the figure boundaries). Below this line, in halo and galactic disk models, 
only non-impulsive models are allowed, where rj and the duration are independendent 
parameters not connected by eq.(2.1.4). The duration could, in this case, be given, e.g., 
by the wind duration ty^ discussed in §2.1, or by an accretion time tacc- For non-impulsive 
models, in Figs. 1 through 9, the duration [ti) values plotted would not be the real 
duration; the "duration" values in Figs. 4-9 can, however, be simply re-interpreted as a 
quantity proportional to an inverse power of rj = Eq/MqC^ , a characteristic of the model, 
while the real duration is an extra model parameter. The cosmological models, however, 
can satisfy the impulsive approximation throughout the region of parameter space here 
considered, so for these the duration of eq.(2.1.4) can be self-consistently determined 
from 77. (Non-impulsive models could also be considered in the cosmological case, at the 
price of introducing a modcl-dcpcndent timescale, but this is not necessary). 

There is a well-defined density dependence of the fiuences, as seen from Figs Id, 
2d and 3d for the piston model (this holds also for the other two models). The reason for 
this is that the deceleration radius (2.1.1) depends on density, and this determines both 
the duration (2.1.4) of the burst and the field strength (2.2.1) at the shock. While the 
break energy behaves monotonically with no, the fiuence does not, since that depends on 
the ratio of synchrotron to IC efficiencies which is a more complicated function of the 
density. The density effect is also seen in Fig. 9, where the lines of — constant and 
Sj/Sx = constant are plotted in the Hq, tL plane. This defines the regions where both 
an "X-ray paucity" constraint and a minimum 7-ray fluence constraint can be satisfied for 
the three models considered, in both the galactic and cosmological cases and for various 
values of the field generation parameter ^ (or A) and the particle acceleration parameter 
K. It is seen that, in general, an external density Ug cm ^ is required to produce 
acceptable fiuences. This would be an argument against having GRB in a galactic halo. 
In the case of GRB in our galactic disk it would predict that they are seen only inside 
spiral arms, where the ISM density is the standard value no ~ 1 cm~^ or greater (but the 
X-ray paucity constraint and burst duration constraint cannot be easily satisfied, unless 
additional assumptions are made, as discussed above). Strong 7 emission could arise if 
the electron spectrum were flatter (with p < 3), so that synchrotron emission from the 
highest-energy electrons was the dominant process (cf. Begelman et al. 1993). In the case 
of cosmological bursts. Figure 9 suggest that the GRB occur in (distant) spiral galaxies or 
other "galactic" environments where the density is of the order or above the standard ISM 
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density of 1 cm~^. However, it is possible that cosmological GRB which have wandered 
out of their host galaxy, or galactic GRB in a halo, could be inside a "local high density 
pocket" produced, e.g. by a precursor wind, and might therefore have fluences comparable 
to those of bursts in a high density ISM. 

The spectral calculations are used also to make predictions of the optical/UV and 
the GeV/TeV fluences, as well as the X-ray fluences expected from bursts with a given 
MeV-band fluence and duration. These fluences can be compared with the limiting sensi- 
tivities of BATSE (e.g. Fishman, 1992) Sy^m ~ 10~^ ergs cm~^ in the 7-ray band, and of 
HETE (Kicker, 1992), S^^^/Il ~ SxlQ-^ ergs cm'^ s'S S^^^/Il ~ 8x10"^ ergs cm'^ s"^ 
and Su,Tn/tL ~ 10~^ ergs cm~^ s~^ in the 7-ray, X-ray and 0/UV bands respectively (lOd 
limits). For a typical burst duration ti, ~ 10 s, the HETE fluence sensitivity limits are 
5^ ~ 3 X 10"'^ ergs cm~^, 5'a; ~ 8 x 10~* ergs cm~^ and Su ~ 10~^ ergs cm~^. 

Signiflcant X-ray fluences are predicted in some of the models. For a given (ob- 
servable) 7-ray fluence, even if one imposes restrictions on the model parameters to satisfy 
an X-ray paucity constraint (e.g., based on current observations, S^/S^ ^ 3 x 10~^), the 
models predict X-ray fluences that ought to be detectable. With current detectors this is 
possible if the GRB happens to be in the (usually rather narrow) X-ray detector field of 
view. The latter, unfortunately, is statistically improbable. For instance, the Rosat field 
of view is 3 square degrees, i.e. ~ 10~^ of the entire sky, so for a burst rate of one per 
day, one would have to wait ~ 10^ days to find a GRB in the field of view. This situation, 
however, should greatly improve with the advent of a new generation of omnidirectional 
X-ray detectors, e.g. HETE (Ricker, et al. , 1993). The predicted X-ray fluences, even 
when satisfying the X-ray paucity constraint, can reach values Sx ^ 10~^ — 10~^ ergs cm~^ 
for the brighter bursts, which is above threshold values for typical detectors, if the source 
is in the fleld of view. 

The optical/UV fluences, again for sources satisfying the canonical X-ray paucity 
ratio, are predicted to be in the range Su ^ 10"^*^ ergs cm~^, and in some cases (e.g. 
the frozen-in cosmological model) as high as 5^ 3 x 10~^ ergs cm~^ for D28 — 1, 
when S'^ ~ 3 X 10~^ ergs cm~^. The nearer (brighter) bursts of ~ lO"'^ ergs cm~^ 
therefore could have 5'^ ~ 3 x 10~^ ergs cm~^, well above HETE's sensitivity. Other 
models, e.g. the piston model satisfying both 7-ray observability and X-ray paucity, would 
predict significantly less UV fiuence, e.g. for 5"^ ~ 10~^ ergs cm~^ it predicts Su ~ 
3 X 10~^ ergs cm~^. The HETE 0/UV band limiting sensitivity can also be expressed in 
U-band magnitudes as my^rn. ~ — 2.51ogio('S'M,m^^M ''^^L^/l-QxlO"^^ ergs cm~^ Hz) ~ 11. 
for a burst of Su/th ~ 10~^ ergs cm~^ s~^. The U-magnitude for an arbitrary burst is 
m„ = 11. — 2.5 log]^o('S'u/10~^ ergs cm~^) -|- 2.51ogiQtL, where tL is the observed burst 
duration in seconds. 

The X-ray and UV fiuences discussed above are appropriate for the case when a 
restriction of the type Sx/ S^ ^ 0.03 is assumed valid, which requires in general a break at 
some energy between the 7-ray and X-ray energies. From the Band, et al. (1993) analysis, 
it appears that many of the bursts observed with the spectroscopy detector have breaks at 
energies ranging from 400 to 100 KeV, and some above that. However, a number of GRB 
do not show a break within the detector energy range (down to ~ 50 keV). This implies 
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that there could be some bursts which have a higher X-ray fiuence than that mentioned 
above, e.g., possibly as much as 5'2;/5'^ ~ 0.1 — 1. If the spectral breaks in some bursts occur 
in the UV, the 0/UV fiuence could then be much larger than for those bursts satisfying 
an X-ray paucity restriction. 

We note that the self-absorption cut-off would, in all our models, suppress radio- 
band emission during the burst, especially in "cosmological" models. There is, however, 
the possibility of a faint radio "afterglow" , delayed by a few days in the cosmological case, 
by which time the surface area of the fireball has become much larger (cf. Paczynski and 
Rhodes 1993). The delay would be shorter for "galactic" models, however, and radio flashes 
that might be observed could then more plausibly be attributed to coherent pulsar-type 
emission arising from violent magnetospheric effects around a neutron star. 

The fluences at all energies depend, of course, on the assumed electron power law 
index and the slopes above and below the break Sa, si,- The extrapolation of our results 
for slope values other than those used here is straightforward, since the break energies and 
the fiuence at these break points do not change. For the fiducial power per decade slope 
above the break assumed here, the GeV fiuences are usually comparable to the MeV fiuence 
(when the MeV band is near or below the main break). The predicted GeV (or higher 
energy) fiuences could also depend on a high energy cutoff of the electrons. (Above photon 
energies ^ 1 TeV photon-photon opacity effects are likely to reduce the fiuences to values 
below those shown in the figures). The EGRET fiuence sensitivity above 1 GeV (Gehrels, 
et al. , 1991) is of the order S oeV ^ 10~^ ergs cm~^ (this is a rough extrapolation from the 
steady source sensitivity for a 5 x 10^ s integration; for a burst, the background would need 
to be evaluated more carefully). The instrument has a ~ 0.6 sr field of view as compared 
to 47r for BATSE. If the power slope is indeed 0, some of the brighter (McV) bursts should 
be detectable by EGRET, and indeed they have been (e.g. Schneid, et al. , 1993, Kwok, 
et al. , 1993). 
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Figure Captions 

Fig. 1: Cosmological model spectra for a luminosity distance Dl = 10^® cm. Eg = 
10^^ ergs, 9 = 10"^ and r] = 10^. The first three arc the cosmological frozen-in, turbulent 
and piston models for an external density uq = 1 and duration tj, = 5 s (top left, top right, 
bottom left). Identifying models by a label (type-distance, logE'o, log6', log no, log^, logn), 
these are FC51-1003, TC51-10-63, PC51-10-32. The fuU lines are the total spectra, and 
the dashed lines are the spectral sub-components discussed in the text. The fourth spec- 
trum (bottom right, PC51-ln-32) is the cosmological piston for three different densities, 
no = 10~^, 1, 10"^, and durations tL = 50., 5., 0.5 s (left to right). A similar density 
dependence also holds for the frozen-in and turbulent models. 

Fig. 2: Halo model spectra for a distance -D = 50 Kpc, 9 — 10~^. The first three, 
for no = 10-^, are the frozen-in (top left, FH45-1-303) turbulent (top right, TH41-1-303), 
and piston models (bottom right, PH41-1-303), for two durations, ti = 10, 0.5 s (dashed 
and full lines). The bottom right figure is the piston halo model PH41-ln03 for densities 
no = 10~^, 1, 10^ and durations tL = 10, 1, 10"^ s (left to right). 

Fig. 3: Galactic disk model spectra for a distance D — 1 Kpc, 9 — 10~^, showing for 
no = 1 the frozen-in (top left, FD43-1003), turbulent (top right, TD39-1003) and piston 
(bottom left, PD39-1003) models at two durations = 5, 0.2 s (dashed, full lines). In the 
bottom right figure the piston disk model PD39-ln03 is shown for three external densities 
no = 10^, 1, 10-2 and durations tL = 2, 0.2, 0.02 s (left to right). 

Fig. 4: Frozen-in cosmological model FC51-1003, for an external density no = 1 
and Dl — 10^^ cm. The model parameters are FC51-1003, giving (type-distance, 
log£^o, log6', logno, log^, logre). The lower left flat portion is not part of the spectrum, 
it represents the HETE threshold at fluence F ~ lO"-*^*^ ergs cm~^. In order of increas- 
ing energies, the breaks arc due to synchrotron and IC in the reverse shock, and reverse 
synchrotron photons IC-scattered in the blast wave. 

Fig. 5: Turbulent cosmological model TC51-10-63 for no = 1, Dl = 10^^ cm. As before, 
the detection threshold is at 10~^^ ergs cm"^, in order of increasing energy the breaks are 
due to synchrotron and IC in the reverse shock, reverse synchrotron photons IC-scattered 
on the blast wave, and blast wave synchrotron photons IC-scattered in the blast (the blast 
synchrotron component is submerged by the others). 
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Fig. 6: Turbulent halo model TH41-1-303 for an external density uq = 10~^ and D = 50 
Kpc. Detection threshold is 10~^*^ ergs cm~^. The main breaks seen arc the reverse 
synchrotron , blast synchrotron and combined IC (barely seen at low durations), with 
increasing energies. 

Fig. 7: Piston Cosmological model PH51-10-32 for an external density no = 1, and 
detection hreshold at 10"-'^° ergs cm~^. The two breaks are due to synchrotron and IC 
photons from the blast. 

Fig. 8: Piston (galactic) disk model PD39-1003, for an external density no — 1, distance 
D = 1 Kpc and detection threshold 10~^*^ ergs cm~^. The main break is due to the blast 
wave synchrotron photons, while the blast wave IC is at very high energies, and relatively 
weaker. 

Fig. 9: The 7-ray band (0.5 MeV) fluence contours and the 7-ray to X-ray band fluence 
ratios (0.5 MeV/1 KeV), in the duration-external density parameter space, for the three 
types of model F, T, P (columns) at the three distance scales C, H, D (rows). The 7 
fluence has three contour levels, = 10~^, 10~^, lO"'* ergs cm~^ (dotted, short-dashed 
and long-dashed lines, unlabeled); the j/X ratio has three contour levels, logS-y/ log Sx = 
l.,1.5, 2. (dot-dash, dot- long dash, and dotted, with label r). In the halo and galactic 
disk cases, models below and to the left of the solid diagonal line have rj > Tmax = 
3.3 X lO^Eli^ Tg (see text), i.e., they don't satisfy the impulsive approximation; in this 
region, the duration is not given by eq.(2.1.4) as plotted, but must be specified as an 
additional free parameter. 
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